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Abstract. We establish a Lipschitz stability estimate for the inverse problem consisting in the determi- 
nation of the coefficient u{t), appearing in a Dirichlet initial-boundary value problem for the parabolic 
equation dtu — AxU + a{t)f{x)u = 0, from Neumann boundary data. We extend this result to the same 
inverse problem when the previous linear parabolic equation in changed to the semi-linear parabolic equation 
dfU — Aa;U = F{t, X, o-(t), u{x, t)). 

Key words : parabolic equation, semi-linear parabolic equation, inverse problem, determination of time- 
depend coefficient, stability estimate. 

AMS subject classifications : 35R30. 

Contents 

1. Introduction 

2. Time-differentiability of potential- type functions 

3. Proof of Theorems 1 and 2 
References 



1. Introduction 

Throughout this paper, we assume that Vlis, bounded domain of K" with n'^ 2. Let T > and set 

Q = r2 X (0, T), r ar2, e = r x (o, r). 

We consider the following initial-boundary value problem 

- A:rit + cr(t)/(a;)u 0, {x,t)eQ, 
u{x,0)^h{x), xen, (1.1) 

u{x,t) ^ g{x,t), (a;,<) e E. 

We introduce the following assumptions : 
(HI) / e C^m, h e C2'°(n), g € C2+"'i+f (E), for some < a < 1, and satisfy the compatibility 
condition 

dtg{x, 0) - A,h{x) + a{{))f{x)h{x) =0, a; e T. 
(112) There exists € F such that 

inf \g{xo,t)f{xo)\ > 0. 
te[o,T] 

Under assumption (HI), It is well known that, for a G C^[0,T], the initial-boundary value problem (|l.ip 
admits a unique solution u = u{cr) S C^"'""'^^^ (Q) (see Theorem 5.2 of |LSU| ). Moreover, given M > 0, 
there exists a constant C > depending only on data (that is Q, T, /, g and h) such that ||cr||t^i,oc.(Q ^ M 
implies 

h(f7)||c.+=,iw2(Q) <a (1.2) 
1 
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In the present paper we are concerned with the inverse problem consisting in the determination of the 
time dependent coefficient a{t) from Neumann boundary data d^u{a) on E, where di, is the derivative in the 
direction of the unit outward normal vector to F. 

We prove the following theorem, where B{M) is the ball of [0, T] centered at and with radius M > 0. 

Theorem 1. Assume that (HI) and (H2) are fulfilled. For i ^ 1,2, let ai G B{M) and Ui = u{ai). Then 
there exists a constant C > 0, depending only on data, such that 

||o-2 - '7i|Il~(o.t) =^ ^ \\dtduU2 - 9t9,yiii|lioc.(s) • (1-3) 

Following |COY| . it is quite natural to extend Theorem [1] when the linear parabolic equation is changed 
to a semi-linear parabolic equation. To this end, introduce the following semi-linear initial-boundary value 
problem : 

dtu — Aa-w = F{x, t, a{t), u{x, t)), (x, t) <E Q, 

u{x,0) = h{x), xeil, (1.4) 

u{x,t) — g{x,t), (a;,i)GS 
and consider the following assumptions 

(H3) h e C^'"(H), g e C^+"'^+'^(E), for some < a < 1, and satisfy the compatibility condition 

dtg{x, 0) - A^,h{x) ^ F(0, X, a(0), h{x)), xeT. 

(H4) F e C^{Tix X Rt X X R„) is such that duF and daF are C\ F and daF are with respect to x 
and u. 

(H5) There exist M > and .to G F such that 

inf \daF{xo,t,a,g{xo,t))\ > 0. 

(H6) There exist two non negative constants c and d such that 

uF{x,t,a(t),u) !^cu'^ + d, te[0,T], x eTi, ueR. 

Under the above mentioned conditions, for any a € C^[0,r], the initial-boundary value problem (jl.4p 
admits a unique solution u ~ u{(j) e c2+a,i+ j-q^ ^ggg Theorem 6.1 in jLSU[ ) and, given M > 0, there exists 
a constant C > depending only on data (that is il, T, F, g and h) such that ||c||vi'1.oo(o,t) ^ implies 

||u(cr)||p2+o.i+o/2(Q) s$ C. (1.5) 

We have the following extension of Theorem [T] 

Theorem 2. Assume that (H3), (H4), (H5) and (H6) are fulfilled. For i ^ 1,2, let at G B{M) and 
Ui = u{ai). Then there exists a constant C > Q, depending only on data, such that 

\W2 - o-i!Iloc(o,t) C\\dtd,jU2 - 9t9i.ui||^oo(s) • (1-6) 
Remark 1. Let us observe that we can generalize the results in Theorems]^ and\^ as follows: 

i) In p.ip . we can replace a{t)f{x) by X]fc=i '^k(t)fk{x), where fk, 1 ^ k ^ p, are known. Assume that 
(HI) is satisfied, with f = fk for each k, where the compatibility condition is changed to 

V 

dtg{x, 0) - A,/i(rE) + '^k{0)fk{x)h{x) =0, x&r. 
fc=i 

Therefore, to each (di, . . . , ap) G C[0, T]^ corresponds a unique solution u = u{ai, . . . , Up) G C^^"'^^"/^((5) 
and max{||(Tfe||vi/i,oo(o.T); 1 ^ fc ^ p} ^ Af implies 

\\u{ai, (Tp)||c2 + o.l + o/2(Q) ^ C, 

for some positive constant C depending only on data. 
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Following the proof of Theorem]^ we prove that, under the following conditions '. there exists • ■ • ? G 
r such that the matrix M(t) ~ {fk{xi)g{xi,t)) is invertible for any t G [0,T], 

max hi - (tI\\l^{o.t) C \\dtdi,U2 - dtdt^uiW , 

if aj. G B{M), 1 ^ k ^ p and j = 1,2. Here C is a constant that can depend only on data and Uj ~ 
u{cr{, . . . ,ct|,), j = 1,2. 

a) We can replace the semi-linear parabolic equation in (|1.4p by a semi-linear integro- differential equa- 
tion. In other words, F can be changed to 

Fi{x,t,cr{t),u{x,t)) + / F2{x, s,cr{t — s),u{x, s))ds. 
Jq 

Under appropriate assumptions on Fi and F2, one can establish that Theorem\^ is still valid in the present 
case. 

a) Both in (|l.ip and (11.411 . the Laplace operator can be replaced by a second order elliptic operator in 
divergence form : 

E = V ■ A{x)V + B{x)-V, 

where A{x) — {aij{x)) is a symmetric matrix with coefficients in C^^"(r2), B{x) ~ {bi{x)) is a vector with 
components in C"(f2) and the following ellipticity condition holds 

A{x)^ ■ ^ ^ X\^\^ , ^eR", xeTi. 

Actually, the normal derivative associated to E is the boundary operator d^j^ ~ ^{^) ' ^(a;)V. 

To our knowledge, there are only few results concerning the determination of a time-dependent coefficient 
in an initial-boundary value problem for a parabolic equation from a single measurement. The determination 
of a source term of the form f{t)xD{x), where xd the characteristic function of the known subdomain D, was 
considered by J. R. Canon and S. P. Esteva. They established in |CE86-1) a logarithmic stability estimate 
in ID case in a half line when the overdetermined data is the trace at the end point. A similar inverse 
problem problem in 3D case was studied by these authors in |CE86-2| . where they obtained a Lipschitz 
stability estimate in weighted spaces of continuous functions. The case of a non local measurement was 
considered by J. R. Canon and Y. Lin in |CL88| and |CL90| . where they proved existence and uniqueness 
for both quasilinear and semi-linear parabolic equations. The determination of a time dependent coefficient 
in an abstract integrodifferential equation was studied by the first author in |Ch91-l) . He proved existence, 
uniqueness and Lipschitz stability estimate, extending earlier results by |Ch91-2| . |LS87| . |LS88| . |P085-1| 
and |P085-2) . In |CY06| . the first author and M. Yamamoto obtained a stability result, in a restricted 
class, for the inverse problem of determining a source term f{x,t), appearing in a Dirichlet initial-boundary 
value problem for the heat equation, from Neumann boundary data. In a recent work, the first author 
and M. Yamamoto |CY11| considered the inverse problem of finding a control parameter p{t) that reach a 
desired temperature h(t) along a curve "/(t) for a parabolic semi-linear equation with homogeneous Neumann 
boundary data and they established existence, uniqueness as well as Lipschitz stability. Using geometric optic 
solutions, the first author |Ch09) proved uniqueness as well as stability for the inverse problem of determining 
a general time dependent coefficient of order zero for parabolic equations from Dirichlet to Neumann map. In 
[EOT] and [E08] . G. Eskin considered the same inverse problem for hyperbolic and the Schrodinger equations 
with time-dependent electric and magnetic potential and he established uniqueness by gauge invariance. 
Recently, R. Salazar |Sa| extended the result of |E07) and obtained a stability result for compactly supported 
coefficients. 

We would like to mention that the determination of space dependent coefficient f{x), in the source 
term a{t)f{x), from Neumann boundary data was already considered by the first author and M. Yamamoto 
|CY06| . But, it seems that our paper is the first work where one treats the determination of a time dependent 
coefficient, appearing in a parabolic initial-boundary value problem, from Neumann boundary data. 
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This paper is organized as follows. In section 2 we come back to the construction of the Neumann 
fundamental solution by [It] and establish time-differentiability of some potential-type functions, necessary 
for proving Theorems 1 and 2. Section 3 is devoted to the proof of Theorems 1 and 2. 

2. TiME-DIFFERENTIABILITY OF POTENTIAL-TYPE FUNCTIONS 

In this section, we establish time-differentiability of some potential-type functions, needed in the proof 
of our stability estimates. In our analysis we follow the construction of the fundamental solution by S. Ito 

m- 

First of all, we recall the definition of fundamental solution associate to the heat equation plus a time- 
dependent coefficient of order zero, in the case of Neumann boundary condition. Consider the initial- 
boundary value problem 

' df U = A^u + q{x, t)u, {x, t) ^ fl X (s, Iq), 

VLm.u{x,t) ~ uo{x), X G ft, (2 1) 

^ d^u{x, t) = 0, (x, t) £V X (s, to). 

Here sq < Iq are fixed, s G (so,io)j ""o s^nd q{x,t) are continuous respectively in f2 and in f2 x [s,to]- Let 
U{x, t] y, s) be a continuous function in the domain sq < s < t < to, x £ fl, y Q. We recall that U is the 
fundamental solution of (|2.ip if for any uq G C{Q), 



u{x, t) 



U{x,t;y,s)uQ{y)dy 



is the solution of (|2.1|) . We refer to [It] for the existence and uniqueness of this fundamental solution. 
We start with time-differentiability of volume potential-type function|3- 

Lemma 1. Fix s G (so,io)- Let f € C{fl x [s,to]) be with respect to x, q £ C^{fl x [s,to]) and define, for 
{x,t) GQ.X {s,to), 

f^{x,t;T)=l U{x,t;y,T)f{y,T)dy, t > t > s. 



(2.2) 



Then, f^ admits a derivative with respect to t and 
dp 



df f 

— {x,t;T) ^ J^U{x,t;y,r){Ay + q{x , t)) f {y , t) dy 

U {x, t; z, T')dtq{z, t')U{z, t' ; y, T)f{y, T)dzdydT' . 

I 

f'^{x,t;T)dT, {x,t) G r X (so,to), 



T JnJn 



Moreover, F given by 

F{x,t) - 

possesses a derivative with respect to t. 



dF 
'dt 



{x,t) = f{x,t) + 



dt 



{x, t; T)dT 



j'^^[x,t;r)dT j\\f{.,- 



dr. 



(2.3) 
(2.4) 



^Recall that ii ip = ip{x,t) is a continuous function then the corresponding volume potential is given by 

tp[x,t)= f f U(x,t;y,T)ip(y,T)AyAT. 
Js Jn 
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Proof. We have only to prove (|2.2p and (|2.4p . because (|2.3p follows immediately from 
Let then uq e C^{^1) and consider the function 



u{x,t) = / U{x,t;y,s)uo{y)dy, x e^, s <t <to. 
We show that u admits a derivative with respect to t and 
dtu{x,t)^dti / U{x,t;y,s)uoiy)dy 



U{x,t;y,s){Ay + qix,s))uQ{y)dy 



U{x, t; z, T)qt{z, t)U (z, r; y, s)uo(?/)dzdydr. 

Is JnJn 

We need to consider first the case ua = ivq G C°°(ri). Set 

w{x,t) = / U{x,t; z, s)wo{y)dy, x£Q,s<t<tQ. 
Jn 

Clearly, w{x,t) is the solution of the following initial-boundary value problem 

' dtw — AxW — q{x, t)w — 0, {x, t) £ Q X (s, to), 
^ lmiw{x,t) = wo{x), X ^ fl, 

^ ^^^(a;, t) = 0, (x, i) e r X (s, <o) 

and itii = dtW satisfies 

' dtWi - AxWi - q{x, t)wi = ~dtqw, (x, t) € Q x (s, to), 
< limwi(x,t) = (Aa: + s))wo(a;), a; £ fi, 

^ 9i.wi(a;,t) = 0, {x,t) & T x {s,to). 

Therefore, p.Sp . with if in place of u, is a consequence of Theorem 9.1 of [Lt] . 

Next, let (wo)n be a sequence in C°°(ri) converging to uq in C^(ri) and v{x,t) given by 



w(a;,t)= / U{x,t]y,s){Ax+q{x,s))uoiy)dy 



Is JnJn 

Consider {wn)n, the sequence of functions, defined by 



U{x, t; z, T)dtq{z, t)U (z, r; y, s)uo(y)dzd?/d7 



w„{x,t) = I U{x,t;z,s)'WQ{y)dy. 

We proved that, for any n G N, 

dtWnix,t) ^ / C/(a;,t;2/,s)(Aj^ -1-9(2;, s))wo(2/)dy 



[/(x, i; z, T)dtq{z, t)U (z, t; y, s)^^ (y)dzd2/dT. 



s JnJn 



From the proof of Theorem 7.1 of [It] . 

|C/(a;,t;y,s)|dy <C7e^(*-^), (x, t) G H x (s, to). 



(2.5) 



(2.6) 



(2.7) 



Therefore, we can pass to the limit, as n goes to infinity, in (|2.6p . We deduce that dtWn converges to v in 
C(ri X [s,to])- But, Wn converges to u in C(ri x [s,to]). Hence u admits a derivative with respect to t and 
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holds true. Finally, we note that 



is deduced 
□ 



dtu = V. That is we proved (j2.5p and consequently 
easily from (|2.7p . 

Next, we consider time-differentiability a single layer potential- type functioij^. 

Lemma 2. Fix s € (so,io)- Let f G C{T x [s,to]) be with respect to t ^ [Sjta] with f{x,s) = 0. Define, 
for e r X is,to), 



f\x,t-T) 



Then 



U{x,t]y,T)f{y,T)dxj{y), t > t > s. 
Fix,t) -- 



is well defined, admits a derivative with respect to t and we have 

dF 



dt 



f\x,t;r)dT 
id we have 

^ C l|5t/llL~(rx(s,to)) ■ 



(2., 



L=-(rx(s,to)) 

Contrary to Lemmafl] for Lemma[2]we cannot use directly the general properties of the fundamental so- 
lutions developed in We need to come back to the construction of the fundamental solution of (|2.ip intro- 
duced by [It]. First, consider the heat equation dtu = A^^u in the half space fii = {x = (xi, . . . , a;„); xi > 0} 
in K" with the boundary condition dx^u = on Fi = {a: = (0, X2, ■ ■ ■ , Xn)', {x2, ■ ■ ■ , Xn) G R"^^}. For any 
y = (2/i,2/2,---,2/n), we define f by y = (-yi, ya, • • ■ , Let 

1 1==!^ 

G{x,t) = j—^e-^ 
(Ant) 2 

denotes the Gaussian kernel and set 

Gi{x,t;y) = G{x - y,t) + G{x ~y,t). 
Then, the fundamental solution Uo{x,t;y, s) of 

dtU = AxU, {x,t) e ni X {s,to), 

VLmu{x,t) = uo{x), xGfii, (2.9) 

^d^u{x,t)=0, {x,t) eTi X {s,to) 

is given by 

Uo{x,t;y,s) = Gi{x,t ~ s;y). 
In order to construct the fundamental solution in the case of an arbitrary domain fi, Ito introduced the 
following local coordinate system around each point z g F. 

Lemma 3. (Lemma 6.1 and its corollary. Chapter 6 of [It]) For every point z E T, there exist a coordinate 
neighborhood Wz of z and a coordinate system {x\, . . . ,a;*) satisfying the following conditions: 

1) the coordinate transformation between the coordinate system (xj, . . . , x*) and the original coordinate 
system in Wz is of class and the partial derivatives of the second order of the transformation functions 
are Holder continuous ; 

2) T D Wz is represented by the equation x^ = and f2 n Wz is represented by x^ > ; 

3) let C be the diffeomorphism from Wz to C{Wz) defined by 

C: W, ^ C{Wz) 



{x*i{x),...,xl{x)). 

function = ip{x, t) if 
ip{x,t) = J J U{x,t;y,T)ip{y,T)da{y)d7 



The single-layer potential corresponding to a continuous function ip = ip{x, t) is given by 
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Then, for any u £ C^(Jl) we have 

From now on, for any z G F, we view coordinate system {x\, . . . , x* ) as a rectangular coordinate system. 
Moreover, using the local coordinate system of Lemma |3l for any y = (2/1,2/27 • ■ • 72/n) S '^(Wz), we define 
y = (— j/i, 2/2, . . . , j/,i) and, without loss of generality, we assume that, for any y € C{Wz), we have y € C(Wz). 
For any interior point z of fi, we fix an arbitrary local coordinate system and a coordinate neighborhood Wz 
contained in fl. For any z £ and 6 > 0, we set W{z, 6) = {x : \x ~ < 6} and Sz > such that, for any 
z g we have W{z, 5z) C Wz- 

Recall the following partition of unity lemma. 

Lemma 4. (Lemma 7.1, Chapter 7 of [It]) There exist a finite subset {zi, . . . , z„i} of fl and a finite sequence 
of functions {wi, . . . , ujm} with the following properties: 

1) supp uji C W{zi,Szi), I = l,...,m, and each to i is of class C'^ with respect to the local coordinates in Wzi ; 

2) {uJi{x)'^; 1 = 1,..., m} forms a partition of unity in Q, ; 

3) d,LUi{0=0, / = l,...,m, CgF. 

Let {zi, . . . , z,„} be the finite subset of f2, introduced in the previous lemma. For any k G {!,..., m}, 
let Ck denotes the diffeomorphism from Wz^ to CkiWz^) defined by 

A: Wz, ^ CkiWz,) 

X h^. {xl{x),...,x'^{x)), 

where (x*,...,a;*) is the local coordinate system of Lemma [3] defined in Wz,. For any k e {l,...,m}, the 
differential operator 

dt- A^- q{x,t) 
becomes, in terms of local coordinate system x* = (a;*, . . . , x* ), 



1 " 

^* / / ,N E ^^l (VaJx*)a^ix*)d^y) - qk{x* ,t) 



in £k{Wz,) X (so,io)- Here qk{x*,t) is Holder continuous on Ck{Wz,) x (so,io) and (a*-' (a;*)) is the con- 
travariant tensor of degree 2 defined by 

[a]^ix*)) ^ {JcJCk\x*))f {JcACj:\x*))) , 

with 

According to the construction of [it] given in Chapter 6 (see pages 42 to 45 of [it]), ^a^-' (a;*)^ is of class 

in £fe(ri n Wzf.) and it is a positive definite symmetric matrix at every point x* € CkiWz,)- We set 
{a'lj{x*)) = (Ofc (2:*))^^ and ak{x*) = det(a|j(x*)). Consider the volume element dx* = ak[x*)dx\ . . . dx* 

on Ck{Wz,) and dec' = y'a{Q,x')dxl . . ■ dec,* on Ck{Wz, n F) with x' = {x^, . . . ,x*,). Note that, by the 
construction of S. Ito [It] (see page 45), for any fc = 1, . . . , m, we have 



a'l{Ck{x)) = a'l[Ck{x)), xe^C^Wz,, for i = j = 1 or j = 2, . . . , n, (2.10) 
al'{Ck{x))=a'i^{Ck{x)) = -al'{U{x)), x&Tlf^Wz,, forj = l,...,n (2.11) 

and 

=af ='5,1, ^^Tr^Wz,, j = l,...,n, (2.12) 

where 5ji denotes the kronecker's symbol. For any k G {1, . . . , m}, let Gk{x, t; y) be defined, in the region 

Dk = {{x,t,y); X, ye £k{Wz,), <t <to~ s}, 
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by 

1 °?j(«)(x.-»i)u,-i>j) 

Gkix,t;y) ^ e ^'-^=1 

Next , define H^^{x,t;y) = Gk{C k{x),t; Ckjy)), for fc e {l,...,m} and e 17 ; Hy^^{x,t]y) = 
Gfc(/:fe(a;),t;£fc(y)) + Gfc(/:fc(x),i;/:fc(2/)), for A: e {1, . . . , m}, G T, x €E W^., and y e W,, ; H,,{x,t;y) = 
if x ^ Wzf. or J/ ^ VFz^ . Consider also H{x, t; y), defined in the region 

D = {{x, t,y); X e n, y S H, < t < to - s}, 

as follows 

m 

H{x,t;y) = '^uji{x)Hz^{x,t;y)uJi{y). 
1=1 

As in Lemma 7.2 of [It], we define successively: 

Jo{x, t; y, s) = {dt - - q{x, t)){H{x, t - s; y)), 

Jk{x,t\y,s) = I I Jo{x,t;z,T)Jk^i{z,T;y,s)dzdT, 
Js Jn 

+ 00 



K{x,t;y,s) = ^Jk{x,t;y,s). 



k=0 

Then, following [it] (see page 53), the fundamental solution of p.ip is given by 

U{x,t; s,y) = H{x,t ~ s;y) + / / H{x,t — t; z)K{z,T;y, s)dzdT. (2-13) 

Js Jn 

We are now able to prove Lemma [2] with the help of representation ()2.13p . the properties of H{x,t;y) 
and K{x, t; y, s). 

Proof of Lemma\^ Without loss of generality, we assume that s = 0. Set 

Fi{x,t)= I / H{x,t- s;y)f{y,s)da{y)ds, 
Jo Jr 

F2{x,t) = / / H{x,t~ T;z)K{z,T;y,s)f{y,s)dzdTda{y)ds. 
Jo JrJs Jn 

According to representation (|2.13p . one needs to show that Fi and F2 admit a derivative with respect to t 
and 

\dtF,{x,t)\ + \dtF2ix,t)\ ^CWdtfh^^r.m^)) (2-14) 
for {x,t) G r X (0,<o)- We start by considering Fi. Applying a simple substitution, we obtain 

F,ix,t)^ f f H{x,s;y)f{y,t-s)da{y)ds. (2.15) 
Jo Jr 

Next, for x G F, there exist Zi,...,/^ C {!,..., m} such that x G suppw; for I G {h, . . . ,lr} and 
X ^ suppw; for I ^ {li, . . . , Ir}. Moreover, since a; G F, we have z/^, . . . , zj^ G F. Then, from the construction 
of H{x, t; y), we obtain 



H{x, s; y)da{y) ^ J^^l^h i^)Hz,^ {x, s; y)uji^ {y)da{y) 

1 



(47rs) 



^ r 1 1- (o,h')(^' -h'Xi' -y' ) 

2E / X;.(0,x')-— ^e-^".-^^ ^ XiM)V^M'W 

f-^lR^-i 47rs 2 



fc=i • 
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with, for / e {1, . . . , m}, xi G Cq (£;(supp w/)) such that x/(a;) ~ ^{x)) and with {x'l, . . . , x'^) = (0, x'), 

{y'l, . . . , y^) — (0, y'). Using the substitution y' ^ z' = ^ ^ , we derive 



H{x,s;y)da{y) 



fe=i 

Therefore, 



(2.16) 



C 

|i7(a;, s;y)|dCT(y) —= 
where ao > is a constant. From this estimate, we deduce that 



\Hix,s-y)fiy,t-s)\da{y)^C- 



and 



H{x,s;y)f{y,t- s)da{y) 



HL°°(rx(o,to)) 

c'f./!lL~(rx(0,to)) 



Thus, Fi admits a derivative with respect to 



/ H{x,t-y)f{y,Q)Acj{y)+ / / H{x,s-y)dtf{y,t-s)da{y)ds 
Jr Jo Jr 

and, since f{y, 0) = for y G F, we obtain 

|atFi(a;,t)KC||aj||^^(P,(o^,^)), (x,t)eFx(0,to). (2.17) 
Let us now consider F2. We want to show that dtF2 exists and the foUowing estimate holds: 

|atF2(x,t)KC||aj||^..(P,(o_,^„, (x,t)eFx(o,to). (2.18) 

For this purpose, using the local coordinate system, it suffices to prove 

\dtF2{C^\0,x'),t)\^Ci\\dtfh^^^rxio,to))^ ((0,x'),t)eA(rnW^,,)x(0,to), Z€{l,...,m}. (2.19) 

From now on we set x = ^(0, a;') with (0, x') E Ci{r W,_^) C {0} x R"-i and we wih show (HHH). First, 
note that 

Jo{z, r; s, y) ^{dr - A. - q{z, t))H{z, s;y) 



L'^,,^H,,{CT\z*),t - s-y)u:i{y) 



1=1 



Y^[L\^^„LOi{Ci\z*))]H,,{CT\z*),r-s-y)u:i{y). 



1=1 



According to the results in Chapter 4 of [It] (pages 26 and 27), using the local coordinate system, we obtain 
L[^^.H.ACi\^l.r-s-Ci\yn)^Y.^''\'^^*)-^\'^y*))-dTil^^^^ 



+ [Bi{z*,y*,d,.) + qi{z*,t)]H,, {Ci\z*), T-s- Ci\y*)), 

where Bi{z* ,y* , dz*) is a differential operator of order ^ 1 in z* with continuous coefficients in z*,y* € 
£;(suppa;;). In view of the results in Chapter 4 of [It] (see pages 26 and 27), combining (|2.10p . (|2.1ip . (|2.12p 
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and (|2.16p . applying the substitution y" = , with z* ~ {zl, z') and y* ~ (0, y'), we obtain 



for < s < T < to and z* E Ci{il O Wz^), where, for j — 0, 1, 2, Pj are polynomials and Jq are continuous 
functions, with respect to t, s G {0,to) and satisfy 

d^'df'J'oiz*,T;y",s;vi) dy" ^Ci, < s < r < io, ^i* > 0, < < to, 



max 

1=0,1,2 Jan^i 

"1+02^1 



for some constant C; > 0. We note that 9„j Jq {{z[, z"),T;y" , s;vi) is not necessarily bounded. Indeed, we 
show 



^ —L, 0< vi <to, j = 0,l,2. 



This representation and the construction of K{z, r; y, s) in Chapter 5 of [It] (see pages 31 to 32 for the 
construction in M" and page 53 for the construction in a bounded domain) lead 



K{Cr\z*),r;y,s)daiy) = J] g, e 



Kj {z*,T;y",s;T - s) „ 
(r - s) 2 



(2.20) 



for < s < T < to and z* G £;(i7 n Wz^), where, for j = 0, 1, 2, Qj are polynomials and Kj are continuous 
functions, with respect to r, s e (0,to) and satisfy 



max_^ / la^^a^^Xj (z*,r;y",s;i;i)|dy" G, < s < r < to, ^i* > 0, < ui < to, 
where Cj > is a constant. Furthermore, using representation (|2.20p . we have, for s < r < t < to, 



H{x,t-T;z) / K{T,z;y,s)f{y,s)da{y)dz 
n Jr 

2 m-. 

uJi{xm,{x,t-T;Cr\z*j)xi{zlQj 

j=Q 1=1 



\/t - S 



Kj {z*,T;y",s;T - s) „ 



(r — s) 2 



dy" 



dz* 



with M!J, {(zjf, . . . ,z*) G M"; zj" > 0}. Then, applying the substitutions z" = ^0= and z[ = -;|^, we 
deduce, in view of the form of the functions Kj, the following 

H{x,t~T;z) / K{z,T;y, s)d<T{y)dz 



E 

j=0 " 



g/(a:^t-T;(z;,0^O,T-s) 



ifj((zl,z"),T;y",s;r-s) „ 
dy 



dz"dz[, s < T < t < to, 



-1 (t-s)5 

(2.21) 

for some continuous functions A'q, A'( and H'^ such that A'q, are C^, with respect to s and t, and the 
following estimates hold: 



\Hl{x\t~T; {z[,z"),T ~ s)\dz" ^Q, < s < t < t < to, 



(2.22) 



max / \d^'df^K'{{z[,z"),T;y",s;vi)\dy"i;iCi, < s < t < to, < i>i < to, (2.23) 

j=0,l I I -I 

ai +a2 ^ 1 
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for some constant C; > 0. Repeating the arguments used for (|2.2ip and applying some results of page 31 of 
[It] , we obtain, for < t < to, 



/ / f \H{x,t~T;z)\ f \K{z,T;y,s)\d<jiy)dzdTds^Ci f f 
Jo Js Jn Jr Jo Js 



1 



^CiV] f\t-.sy-ids^Ci. 



drds 



(2.24) 



This estimate and Fubini's theorem imply 



F2ix,t) 



mH{x,t-T;z) / K{z,T;y,s)f{y,s)da{y)dzdTds. 
.1 Jr 



Then, in view of representation (|2.2ip . for all < t < to, 

H'i{x',t-T; {z[,z"),T-s) 

j=o " *+ 



F2{x,t) 



t rt 1 

E 



K'^{{z[,z"),T-y\s;T-s) ^ „ ,; 



(r — s) 2 



-h{x ,s;y ,z )dy' 



dz"dzJ^dTds, 



where fi{x' , s;y" , z") = f (^Ci "'^(0,a;' — {\/t — s)z" — {^Jt — s)y"), s) . Making the substitution t' = t — r, we 
obtain 



F2it,x) 



t rt-s 1 

E 



H'i{x',T';{z[,z"),t- s-t') 



K'{{z[,z"),t-T'-y",s-t-s-T') 



h{x',s;y",z")dy" 



dz"dzJdT'ds. 



Then, the substitution s' = t — s yields 



t ps' 1 



F2ix,t) 



E 



Hlix',T';iz[,z"),s'-T') 



K'^iiz[,z"),t-T';y",t-s';s'~T') , 



fi{x', t - s'; y\ z")d?/"d0;dz"dT'ds' 



But, for < t' < s' < t < to, estimates ([2:221) . ([2:231) and f{y, 0) = 0, y G T, imply 

H[{x',t',x'-{z[,z"),s' -t') f K'^{{z[,z"),t ^T';y",t-.s';s' ~t') 

{s' -T')i 

°°(rx(o,to)) 1 



E 



(2.25) 



xh{x',t-s';y",z")dy"dz"dz[ 



j=0 



{s' ~T')i 



and 



j=o - 



H[ix',T';{z[,z"),s' -t') r K'^{{z[,z"),t-T';y",t-s';s' -t') 
Vt' Jr"-i (s' - T')i 

^, t I I . I II ii\A iiA I A ii\ ^ n ll"5t/llL°°(rx(0,to)) 1 
X /i(x ,t- s ;j/ ,z )dy dz^dz J ^ Q 2^ 



(2.26) 



(s' - T')i 
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From estimates (|2.25p . (|2.26p and f{y,0) —0,yGT,we conclude that F2 admits a derivative with respect 
to t and 



t ps' . 1 
^0 



4j: 



3=0 



H[{x',t';{z[,z"),s' -t') f K'jiiz[,z"),t-T';y'\t-s';s'~T') 

(s'-r')5 



X fi{x' ,t ~ s';y"; z")dy"dz[dz"'jdT'ds' . 



Moreover, (P?^ and (P?^ imply (^1^ and ((TT5)) . Finally, we obtain ([^1^ from (PT7)l and ^J^. This 
completes the proof. □ 



3. Proof of Theorems 1 and 2 

Proof of Theorem [7J Let u = ui — U2 and cr — (J2 — Ui. Then u is the solution of the following 
initial-boundary value problem 

' dtu ~ A^u + (T2{t)f{x)u = a{t)f{x)ui{x,t), {x,t) e Q, 

< u{x,Q) =0, a; e 17, (3.1) 

u(x, t) = 0, (x, t) e S. 

Let U{x,t;y, s) be the fundamental solution of p.ip with q{x,t) = —(J2{t)f{x). Applying Theorem 9.1 
of Ht], we obtain 



u{x, t) 



Jn 



U{x,t;y,s)a{s)f{y)ui{y,s)dyds+ / / U{x,t;y, s)di,u{y, s)da{y)ds. (3.2) 

'0 Jr 



Now, since u{x,t) = 0, (i,.T) G E and x G F, 



[/(x, t; y, s)a{s)f{y)ui{y, s)dyds 







[/ (x, ?/, s)dyu{y, s)da{y)ds. 



(3.3) 







In view of differentiability properties in Lemma [T] and [21 we can take the ^-derivative of both sides of 
identity We find 



f{x)g{x,t)a{t) 



dt[ / U{x,t\y,s)(T{s)f{y)ui{y,s)dy]ds 
\Jn 



dt 



Jr 



U{x,t;y,s){di,U2{y,s) - a,.Mi(y, s))dCT(y)di 



and, for x = xq, condition (H2) implies 



a{t)==h{t)j dtyj U{xo,t;y,s)a{s)f{y)ui{y,s)dy]ds 



+ h{t)dt{ / / U{xQ,t;y,s)di,u{y,s)da{y)di 



Jr 



(3.4) 



where h{t) = -l/{g{t,xa)f{xo)). 

Since u{x,0) ~ 0, x E ft, we have d,ju{x,0) = 0, x G F. Thus, the estimates in Lemma [T] and [2] lead 



dt[ / U{xo,t;y,s)(T{s)f{y)ui{y,s)dy]ds 
\Jo 



dti / U{x,t;y,s)d^u{y,s)da{y)ds 



Jr 



< C / |a(s)|ds, 
Jo 



(3.5) 
(3.6) 
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Therefore, representation p.4p and estimates (|1.5p . p.5p . (|3.6p imply 

|a(t)| rc|a(s)|ds + c 119*9.^^11^^(5.). 
Jo 

Here and henceforth, C > is a generic constant depending only on data. Hence, Gronwall's lemma yields 

\ait)\ s; C 119*9,^11^^(5,) e^* < Ce^^ \\dtd,u\\ , t € (0, T). 
Then (|1.3p follows and the proof is complete. □ 

Proof of Theorem 2. Set u ~ ui — U2 = u{ai) — u{a2). Then, according to (H4) and (H6), u is the 
solution of the following initial-boundary value problem 

' dtu - AxU - q{x, t)u = F {t,x,cri{t),U2{x,t)) - F {t,x,cr2{t),U2ix, t)) , {x,t) € Q, 

< u{x,0) = 0, xen, (3.7) 

u{t, x) = 0, (i, x) € S, 

with 

r-l 



q{x,t) ^ / 9„F [i, x, (Ti(i), 1*2(0;, i) + r(ui(a;, t) — U2(a;, t))] dr. (3.8) 
Jo 

Note that assumptions (H4) and (H6) imply that q e C^{Q). 
On the other hand, in view of (H4), 

F {t, X, ai{t),U2{x, t)) - F (i, X, G2{t), U2{x, t)) = (CTi(i) - a2{t))G{x, t), 
with ^ 

G{x,t)= / 9^F(t,.T,a2(i) + s(ai(i)-(T2(t)),U2(x,t))ds. 

Using this representation, we deduce that u is the solution of 

dtu - A^u - q{x,t)u = i(Ji(t) - a2it))G{x,t), {x,t) e Q, 

u(a;,0) = 0, xen, (3.9) 

u{x,t)^0, (t,x)eT,. 

Let us remark that (H4) and (H6) imply that G G C'^'^{Q). Let U{t,x; s,y) be the fundamental solution of 
(|2.1[) with q{x,t) defined by p.Sp . Then, according to Theorem 9.1 of |It], for (T{t) ~ cri(t) — a'2(i), we have 
the representation 

u{x,t) = / / [/(a;,i; J/, s)cr(s)G(y, s)dyds + / / U{x,t;y,s)d^u{y,s)da{y)ds. 
Ja Jn Jo Jr 

Since u{x,t) = 0, {t,x) € S, we obtain 

U{xo,t;y,s)a{s)G{y,s)dyds = - / U{xo,t;y,s)dt,u{y,s)da{y)ds, (3.10) 
"'o "'0 Jr 

with xo defined in assumption (H5). Combining Lemma 1 and Lemma 2 with some arguments used in the 
proof of Theorem 1, we prove that /i and /2 defined respectively by 

hit) ^ / / Uixo,t;y,s)a{s)G{y,s)dyds, 
Jo Jn 

f2{t) = j J^U{xo,t;y,s)d^u{y,s)da{y)ds, 
admit a derivative with respect to t and 

f[{t)^a{t)G{xo,t)+ f dt( f U{xo,t;y,s)ais)Giy,s)dy]ds, 
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f2(.t)=J dt\^J U{xo,t;y,s)d^u{y,s)d(T{y) ] ds, 



dti / U{xo,t;y,s)a{s)G{y,s)dy ] ds 



|a(s)||lG(.,s)|lc.(o)ds^C / W{s)\ds 



s$ c||9ta^u||^oo(s) 



and 

I dt(^j U{xo,t;y,s)d,yu{y,s)d(j{y)^ ds 

Here and in the sequel C > is a generic constant that can depend only on data. 
Taking the t-derivative of both sides of identity p.lOp . we obtain 

a{t)G{xo,t) = ~ [ dt( [ U{xo,t;y,s)<j{s)G{y,s)dy^ ds 

dt(^J U {xq, t; y, s)d^u{y, s)da{y) ] ds. 
Let us observe that (H5) and niax(||cri||^ , ||cr2||o2) ^ M imply 

\G{xo,t)\= [ \d,F{t,xo,(J2{t) + s(<Ji{t)-<j2{t)),g{xo,t))\ds 
Jo 

^ inf \daF{t,xo,cr,gixo,t))\ > 0. 

tG[0,T],crG[-A/,Af] 



ait)^H{t) dti Uixo,t;y,s)a{s)Giy,s)dy]ds 



Then, 



+ H{t) J dtyj U{xo,t;y,s)d^u{y,s)daiy)] ds, 

where H{t) = -l/G{xo,t). Hence, (|3TT|) . ((3J2l) and (|3T3l) imply 

\<j{t)\^ f C\<j{s)\ds + C\\dtd.u\\^^,^y 
Jo 

We complete the proof of Theorem [5] by applying Gronwall's lemma. 
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